Renormalons as a Bridge between Perturbative and 
Nonperturbative Physics ft. 



V.I. Zakharov 

Max-Planck-Institut fuer Physik 
Foeringer Ring 6, 80805 Muenchen, Germany 

and 

The Randall Laboratory of Physics 
University of Michigan 
Ann Arbor, MI 48109-1120 

Abstract 

In two lectures, we overview the renormalon and renormalon-related 
techniques and their phenomenological applications. We begin with a single 
renormalon chain which is a well defined and systematic way to specify the 
character of corrections in inverse powers of the total energy to observables 
directly in Minkowski space. Renormalons demonstrate also presence of 
nonperturbative contributions. We proceed then to multirenormalon chains 
and argue that they are in fact not suppressed compared to a single chain. 
On one hand, this phenomenon might be a mechanism for enhancement of 
power corrections. On the other hand, the derivation of relations between 
power corrections to various observables becomes a formidable task and asks 
for introduction of models. In the concluding, third part we consider dy- 
namical models for nonperturbative effects, both in infrared and ultraviolet 
regions, inspired by renormalons. 
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1 Introduction 



Renormalons have various facets and the understanding of the renormalon physics went 
through various phases. The very existence of renormalons was realized within a pure field- 
theoretical framework [[[]] and by renormalons one understands usually a simple perturbative 
type of graphs with a chain of many insertions of vacuum bubbles into a bosonic line. The 
renormalons were used first to discuss matters of principle in field theory. In particular, 
the renormalon chain results in n! divergence of the expansion coefficients a n where n is the 
order of perturbation theory (for a review and further references see, e.g., [§, ||). Such a 
behaviour of a n implies in turn that the perturbative expansions are divergent and cannot 
be a full answer. Thus, nonperturbative contributions are called for. 

Renormalon-based phenomenology emerged much later |fj, |3|, |], || and reached its peak re- 
cently (for reviews see, e.g., 0). Essentially, renormalons allowed to specify the character of 
power-like corrections in various observables directly in Minkowski space. Namely, if a phys- 
ical process is characterized by a large mass scale Q (like total energy in e + e~-annihilation 
to hadrons), then the power corrections assume a generic form [Kqcd/QY \ with k integer 
and depending on the variable considered. The power corrections are encoding nonpertur- 
bative contributions. Renormalons are effective to find the value of k in each case. Rather 
loosely, we can call it renormalon kinematics. This kinematical part is acquiring a status of 
something reliable and simple. The applications, however, are limited in nature, as of any 
kinematics. 

A kind of renormalon dynamics shows up first once one tries to pin down more precisely 
what renormalons actually mean. It turns out that in fact the class of renormalon graphs 
is poorly defined 0. Namely, many renormalon chains are the same important as a single 
chain. This observation calls for reconsideration of the results obtained in a single-renormalon 
chain approximation. Generically, the power k is not changed but the relations between non- 
vanishing contributions are difficult to derive in a consistent way. In particular, the model 
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with the running coupling frozen when applied to power corrections |§ appears to be too 
restrictive since it does not allow for many renormalon chains. 

Thus, at the current stage, we are invited to invent dynamical models for non-perturbative 
effects which incorporate or are inspired by inclusion of many renormalon chains. At the 
moment, it is rather a fascinating challenge than a systematic program. 

Now, a few words on the outline of this text. The lectures at the School were intended to 
be pedagogical in nature. However, it is rather hopeless to expose the renormalon technique 
in sufficient detail within the format of this text. Thus, we address the reader to the reviews 
H and references therein for a systematic introduction to renormalons. Here, we concentrate 
mostly on the problems encountered in development of the renormalon physics and on the 
hints for the nonperturbtive dynamics of QCD provided by renormalons. 

2 Single renormalon chain. 

Infrared-sensitive contributions. 

We start with a brief overview of the renormalon routine. The general problem which 
the renormalons are invoked to solve can be formulated as follows. Imagine that we consider 
a physical process with large characteristic mass scale Q, Q ^> Aqcd, so that the running 
coupling a s (Q 2 ) is small. Then if we evaluate an infrared safe variable, it is given as a per- 
turbative expansion in a s (Q 2 ). This is the standard perturbative QCD. The corresponding 
Feynman graphs, while dominated by momenta k ~ Q include also integration over kine- 
matical regions which are characterized by small momenta k ~ Aqcd- Then the effective 
coupling is in fact large and there are no reliable ways to evaluate the contribution of these, 
infrared-sensitive regions. It is apparently contaminated by nonperturbative effects as well. 
The question is, how big is this non-controllable contribution. It is rather obvious that 
the price for focusing on non-characteristic momenta, k -C Q, is a power-like suppression, 
{Kqcd/QY ■ Indeed the corresponding phase space is small. This can just be considered 
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a definition of what we mean by saying that the characteristic momenta are of order Q, 
not Aqcd- One can then try simple-minded estimates of the contributions coming from 
the infrared-sensitive "corners" of the phase space. Renormalons promote such estimates 
to a status of science. The power corrections, in turn, can be crucial at the pre-asymptotic 
energies, or for the precision tests of QCD. 
Example of the thrust. 

For the sake of definiteness we will consider linear corrections (which are leading ones) 
in shape variables [[5], |^, |TD[. We follow Ref. |J where further details of the derivation 



outlined below can be found. 

Consider the thrust variable which is defined as: 

rr Si \Pi,n\ / 1 

T = max n (1) 

£» \Pi\ 

where Pi >n is the component of ith particle momentum, pi in the reference direction n which 
is chosen to maximize the right-hand side of above. Thus, for a twojet event, T = 1 and 
for a completely symmetric 3-jet event T = 2/3. We will consider the simplest final state 
which gives rise to a non-trivial prediction for the thrust, namely, quark- ant iquark pair and 
a gluon, qqg. Then for soft gluons with uj ~ Aq C d we estimate deviations from T = 1 as 

h rp\ f d ^uj A Q cd . . 

(l-T) soft ~ / a s — — ~ a s . (2) 
Jsoft uj Q Q 

Here, the proportionality to a s ■ (du/u) is the standard factor due to a soft gluon emission, 
while the factor uj/Q is due to definition of the thrust since J2i \p%\ ~ Q- We cannot specify 
within this simple estimate the argument of the strong coupling constant a s . It is natural 
to assume, however, that for soft gluons a s is of order one. 

While the very existence of the 1/Q corrections can be demonstrated in a simple way, the 
evaluation of the renormalon chain, with all the coefficients fixed, is much more complicated. 
It goes, however, through a well defined sequence of steps. First, one considers one-loop 
perturbative calculation. 



Note that the virtual gluon correction drops off from (1 — T) because for a two particle 
final state T = 1 and we consider, therefore, real gluon emission, e + e~ — > qqg. For the final 
state considered the thrust is known to be 

T = max(xi,X2,X3) (3) 

where Xi = 2Ei/Q and E% are energies in the CM frame. For massless quarks and gluons 
one finds the range of the variables as: 

< X2 < 1; 1 — Xi < xi < 1; x± + x<i + x% = 2.. (4) 

After a cumbersome algebra one can find that the 1/Q contribution to the thrust comes 
from the following term: 

^->r d ii:^-^r' 2 

where x = k\j Q 2 and precisely this contribution gives the logarithmic term [ln(l—T)]/(l—T) 
in the thrust distribution. Moreover, the 1/Q term is associated with small x, as anticipated. 
For small x the limits of integration in Eq. (||) can be represented as an expansion: 

T x = l-v^+..., T 2 = l-Ax + ... , (6) 

and we need to keep only the first terms in the expansion. 

This calculation is still only a preliminary step to insert the renormalon chain. The 
crucial point to proceed further is that a s is in fact known to run as a s {k 2 1 ). Therefore, it is 
useful to rewrite the expression for the 1/Q correction as: 

Infrared renormalon in the thrust. 

Having derived expression (|7|), one can identify the renormalon in various ways. Note 
first that the integral over k±_ in Eq. (0) is obviously linearly divergent as a function of the 
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upper bound Q so that, upon dividing by Q, the expression (|7|) is of order a s {Q 2 ). But this 
contribution is absorbed already into the standard first-order radiative corrections to the 
thrust. Now we are interested in parameterizing the IR sensitive contribution to (^). 
To this end we use the Borel representation for the running coupling: 



f°° ( ( k\ 

a s {k 2 1 ) = J daexp ( — a ( b ln-^- 



QCD / / 

2„ 



where we assumed for simplicity one-term (3 — function, 60 = (11 — ^nf)/4ir. As the next 
step, substitute (j|) into (|7|) and integrate first over kj_. The integral over kj_ produces then 
a pole in the cx-plane which is due to the renormalon: 

(iy**(*i)*x) = (9) 

V^s " / renormalon ^ pole Q 

To parameterize the IR contribution we are free to define the integral over pole as, say, its 
principal value. In this way we come to 

\h I a s{k\)dkA = const ^ Q Q D , (10) 

V ^* " / renormalon ^ 

where the value of the const depends on the definition of the integral over the pole and the 
procedure described can be called Landau-pole parameterization of the power-like corrections. 
Another popular way to deal with the renormalons is to re-express a s {k 2 ± ) in terms of 

K i - m.( wcyAl) = ?M<ra>)"o™ (ii) 

where we used the same approximation of a single-term /3-function. The reason to re-express 
everything in terms of a s (Q 2 ) is that eventually we are studying the perturbative expansion 
in a s (Q 2 ) for the thrust. Upon substituting (|TT|) into (^) we come to the integrals of the 
type 

[ Q ln n (Q 2 /kl)d^ « Tn\ . (12) 
Jo Q 
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Note that at large n the integrand has a sharp peak at k\ ~ e 2n Q 2 . 
Thus, renormalons produce a perturbative expansion of the type 

(l-T) ~ E «n(a s (Q 2 )) n ~ E (2b ) n n\(a s (Q 2 )r (13) 

Zarge n ;arge n 

which is clearly divergent and can at best be an asymptotical expansion. By using the 
standard techniques, it is straightforward to see that the expansion ( |13D cannot approximate 
a physical quantity to an accuracy better than of order (Aqcd/Q)- Expansion of ( IT]) in 



a s (Q 2 ) looks so as if we have a chain of bubble insertions into the gluon line. Thus, the 
procedure can be called renormalon- chain parameterization Q of the power-like corrections. 
Since the contribution flT^) to a n comes from effective k 2 ^ <C Q 2 it is called infrared (IR) 
renormalon. 

Kinematical nature of power corrections. 

The beauty of renormalons is that, being a pure perturbative construct, they tell some- 
thing about nonperturbative contributions as well. Indeed, the divergence of the perturbative 
expansion at large n, see Eq. (fL3|), implies that a unique result for the physical quantity can 
be obtained only with account of nonperturbative contributions. In this way we learn that 
the thrust receives nonperturbtaive contributions of order Aqcd/Q- Renormalons fix in a 
way a minimal set of nonperturbative terms. 

On the other hand, the renormalon chain is a simple perturbative graph and as such 
is not necessarily a unique device to probe infrared regions. Thus, we can expect that the 
fixation of the value of k characterizing the power correction is a matter of kinematics and 
gauge invariance. Indeed, the simple estimate of the infrared-sensitive contribution to the 
thrust given above did reproduce the main result that the correction is linear in 1/Q. 

2 The example given can actually be somewhat misleading in the sense that the graphs with the vacuum 
bubble insertions depend rather on the virtuality of the gluon, k 2 than on k\ . The k\ as the argument of 
the running coupling emerges in the case considered after further algebra involving the emission probability 
(see, e.g., the second paper in Ref. A more literal use of the renormalon chain can be found in the next 
section. 
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More systematically, one can utilize a finite gluon mass as a probe of infrared sensitivity 
||. Namely, the general rule [II] is that terms non-analytical in the gluon mass squared, A 2 
are infrared sensitive. Moreover, the correspondence between the non-analytical terms and 
the power corrections in real QCD reads as: 



7 \2 'v/A 2 A 2 2 A 2 \/A A 4 2 , . 

aoa s lnA + a\a s — — — \- a,2a s — LnX + asa s — — h a±a s — LnX + ... (14) 



Q S Q 2 Q 3 S Q 4 
A A 2 A 3 A 4 
OqITIAqcd + Oi— — h °2— H 3 — — h 4 — 1" ... 



where the coefficients Oj are calculable in terms of one-loop Feynman graphs while the cor- 
responding coefficients bi incorporate the effect of large-scale fluctuations, both perturbative 
and nonperturbative, of the real QCD. Note that a , b Q = for an infrared safe quantity, as 
is always assumed to be the case. 

In particular, evaluation of the 1/Q correction to the thrust within the model with finite 
A 2 reveals a typical source of non-analytical in A 2 terms, that is the boundaries on the 
the phase space. To find a s ■ (X/Q) corrections to the thrust one proceed along the same 
lines as above. Namely, the value of the thrust for a particular momentum configuration is 
determined by values of x,, = in the CM frame. A crucial point is that for a nonzero 
gluon mass A the range of the variables contains X/Q terms: 

< x 2 < 1-A (15) 
A l-xs-A/Q 

1-^2-7: < Xt < " . 

Q l-x 2 

Furthermore, one computes the thrust distribution (l/ao)da/dT by substituting the phase 
space region according to 2^2,3 being the largest fraction and finding (l/ao)da/dT in each 
region. Adding these, gives the total thrust distribution from which < 1 — T > is easily 
obtained. The result is 0: 

(1-T) 1/Q = 4CV^- (16) 



The procedure described can be called finite- gluon-mass parameterization of the power-like 
corrections. 

It is worth emphasizing that the contribution to the 1/Q term comes from the region of 
energetic quarks and antiquarks only and, in these regions, from the phase space constraints 
(16). It is therefore seen to be soft-gluon dominated, i.e. coming from gluon with energies 
of order A. Moreover, one can establish a direct correspondence between the 1/Q correction 
to the thrust evaluated via the renormalon chain (considered above) and via introduction of 
finite A. Namely, the identification |5j 

a s (X)X -> n Q a s (k 2 ± )dk± (17) 
Jo 

brings the X/Q corrections to the form found within the renormalon technique. 

This example, as well as many others, seems to indicate strongly that the fixation of 
the index k of the power corrections in terms of renormalons is in fact based exclusively on 
the consideration of phase space and of gauge-invariance constraints. Note, however, that 
while the calculation with a finite gluon mass is a straightforward one-loop calculation, the 
clarification of the argument of the running coupling asks in fact for evaluation of two loops 
and may be a subtle matter. Imagine, for example, that the argument of a s would be the 
invariant mass of the quark- gluon system a s (m 2 qg ) rather than a s (k 2 J. Then for soft gluons 
m 2 qg ~ Aq CD only for u g ~ Aq CD /Q <C Aqcd, not u g ~ Aqcd as discussed above. Then, 
there would be no linear corrections at all. 

Open questions. 

There exist at the moment open problems which arise already on this, kinematical level 
of the renormalon calculus. In a way, they are reflections of the same question, whether 
different infrared probes give the same answer. We will spell out in some detail two open 
questions: 

(a) In a single renormalon approximation, power corrections can be studied in case of deep 



inelastic scattering as well |12| . In view of the richness of the data, these predictions may 



provide an important test of the theory (for a review and further references see the third 
paper in Ref. ||). However, there exists an unresolved controversy over the predictions 
themselves In more detail, one considers, as usual, the amplitude for 7* + q — > q ' + g 



where 7* is the virtual photon with 4-momentum q and q is the parton of momentum p but 
keeps the gluon mass A finite. Then the leading non-analytical in A 2 contribution to, say, 
the longitudinal structure function Fl takes the form |13| : 



(FLh , = c^y-^f, (18) 

where z = Q 2 /2(p-q), q 2 = —Q 2 . If one replaces the X 2 ln\ 2 term by {Aq CD /Q 2 ) according to 
the rules discussed in the preceding subsection, then the result coincides with the evaluation 
of the renormalon chain performed earlier [JT2]. The most interesting point about ( |T8D is 
that the z dependence is fixed and this z dependence of the power correction can be checked 
experimentally. 

On the other hand, if one uses the observation that the coupling runs as a s {k 2 1 ) and 
applies the same technique as above then the resulting z-dependence is different: 

{Fl)iIq2 ~ 2^—z-Q 2 - 1 (19) 
where the explicit expression for the integral / can be found in the original paper [pf 



While the both derivations are straightforward, no understanding of this discrepancy was 
reached so far. 

(b) In some cases the leading 1/Q corrections cancel. The best known and first discovered 
14| example of this kind is the Drell-Yan (DY) process. The difference between, say, thrust 



and the DY cross section can be traced [0 to the difference in resolutions needed to measure 
the corresponding observables. Namely, in case of the thrust the momenta of individual 
particles are assumed to be measured with the precision of order Aqcd while in case of the 
DY process the resolution does not involve Aqcd- The proof of this assertion [HJ utilizes 
the Kinoshita-Lee-Nauenberg and Low theorems. Now, an unsolved problem [|l(| is that 
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inclusion of higher loops brings mixture of soft and collinear divergences . As a result the 
Low theorem holds naively only for frequencies u < m 2 /E where m q is the quark mass and 
E is its energy. In the abelian case it can be extended to u ~ m q [[T7J]. In the nonabelian 
case such an extension has not been proven so far. Thus, in case of the Drell-Yan cross 
section there is no proof that the index k does not change from k = 2 to k = 1 once two-loop 
corrections are included. 
Summary to part I. 

We can summarize the presentation above by saying that there exist well defined and 
systematic techniques to parameterize infrared sensitive contributions to the infrared safe 
quantities in terms of the power corrections (Aqcd/QY '■ Despite some open questions the 
techniques might well find their way to textbooks soon. In the next section we will argue, 
however, that the renormalons are not exhausted by a single chain. 

3 Multi-Renormalon Chains. 

In perturbative QCD observables are calculable as a series in a s (Q 2 ), f = Y J n a n{.ct s {Q 2 )) n . 
As far as a s (Q 2 ) is small enough, each next term in this expansion is small compared to 
the previous one (we should qualify this statement, in fact, for n to be not too large). The 
power-like corrections revealed by renormalons can be viewed as representing kinematical 
regions near the boundaries of the phase space. One might think about them as of a small 
part of the corresponding perturbative contribution. However, the basic observation on the 
multi-renormalons is that the power-like corrections associated with higher loops of the per- 
turbation theory are in fact not suppressed at all. One may say that the whole perturbative 
expansion collapses once projected onto the power-like corrections. This phenomenon was 
observed first on the example of ultraviolet renormalon M and reiterates itself in other cases 



as well JT8|, [L9|, |20| (for related discussions see also pl|). In the subsequent subsections 



we will consider examples of the heavy-quark potential at short distances pOJ, ultraviolet 
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renormalons , dispersive approach to the running coupling [ 19 ] . All these examples are rel 
evant also to discussion of dynamical models in the next section. 
A simple estimate. 

The origin of the phenomenon of the "collapse" of the renormalon chains can be under- 
stood in terms of the simple estimate of the soft particle contributions to the thrust given at 
the beginning of the previous section. Namely, let us consider now two soft gluons instead 
of a single gluon. Neglecting the nonabelian nature of the gluons, we may then write: 



(1-?W ~ / « s - / a s — ^--^ ~ (20) 

J soft u>\ UJ\ J soft u)2 ^6 

i.e., the correction is of the same order as for a single soft gluon (see Eq. @). Two effects are 
combined to escape possible suppressions. First, the effective coupling constant is of order 
unit since we are dealing with soft emission. Second, the origin of Kqcd/Q correction is pure 
kinematical. Namely, the deviations from T = 1 because of the emission of a soft particle 
is of order u/Q. This is the price for considering an unnatural-scale momenta, u ~ ^qcd 
(and not uj ~ Q). It is crucial therefore that this price is paid only once, independent of the 
number of soft particles. 

Heavy-quark potential at short distances. 

A weak point of the estimate in the preceding subsection was the assumption that the 
coupling is of order unit. One can circumvent this difficulty, as we will argue here on the 
example of the heavy-quark potential at short distances |20| . 



As the distance r between the quark and antiquark tends to zero, r <^ Aq^ d , the potential 
becomes Coulomb-like: 

\imV(r) = -^5> n ^(r) + ]>> n A^V n - (21) 

' T n n 

where a n , v n are expansion coefficients (a\ = 1) and we reserved for power-like corrections 
as well. 
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To pick up the renormalon contribution, it is convenient to start from the momentum 
representation. To first order in the running coupling: 

-.r, x ^ f ^ 3 k fA 9s s exp(ik • r) . . 

V (r) = -C F J ^(H^ 2 )) P \ 2 (22) 



representation, a s (k ) = / dcr f-r J , to get 



where (only in this subsection) k 2 = k 2 . As the next step we are using again the Borel 

„ M 2C F /•« , f°° „ sin(kr) f k \' 2ab ° . . 

F r = / dor / rffc 1 ; . 23 

n Jo Jo kr \AqcdJ 



Moreover, since we are interested in small-r behaviour we expand the sine in kr and get 
for the corrections to the Coulomb-like potential: 

SV(r) = v A Q cD + v 2 A 3 QCD r 2 + ... . (24) 

It is noteworthy that odd powers of r are absent for this single renormalon chain. We shall 
come back to this point in the next section. 

Our central point now, however, is the contribution of many renormalon chains. This can 
be extracted in a remarkably simple way once perturbative expressions for V(k 2 ) are available 
in higher loop approximation. In fact first three terms in the purturbative expansion were 
calculated explicitly 



2 = _ 4C F 7ra MS / 1 + %s (2 5g3 _ Q 27g ) + ^MS( 39 , 650 _ A1Aln + 0.077^) j 

(25) 



Moreover, in the approximation of the one-loop /3-function, 



a l(k 2 ) = ^r^QCD-n^a s {k 2 ). (26) 

20 (LAqcd 



and by differentiating n times Eq. (24 ) with respect to Kqcd we immediately find con- 
tribution of (n + 1) renormalon chains, associated with the a™ +1 (k 2 ) in the perturbative 
expansion. 
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Thus, we have a unique possibility to compare the contributions of one-, two- and three- 
renormalon chains: 

5V{r) « (const)A 3 QCD r 2 (l + l.l + 6.0 + ..). (27) 

. However, there is no convergence of the series in sight in that case either. 

To summarize, we have shown that the power corrections associated with many (infrared) 
renormalon chains are not suppressed compared to the power corrections introduced via a 
single renormalon chain. 

Coupling freezing vs renormalons. 

The picture of infinitely many renormalon chains, all of which are equally important, is 
difficult to be transformed to a well defined model. The simplest assumption is that actually 
a single chain reproduces correctly the nonperturbative dynamics of QCD. There are well 
known models to this effect. We have in mind, first of all, the model with freezing of the 
running coupling in infrared at a low value. If the effective coupling is always small, then two 
soft particle emission (see the discussion at the beginning of the section) is suppressed and 
the major uncertainty in theoretical predictions disappears. One may relate then various 
1/Q corrections in much more straightforward way 0. In abelian case, one can use the model 
with a finite gluon mass (see section 2) in higher-loop approximation. Since the emission of 
soft abelian gluons is independent this model amounts again to the dominance of a single 
renormalon chain. 

Thus, one may say that these models contradict our relations between, say, one and two 
renormalon chains. It is easy to understand the reason for this discrepancy. Namely, the 
renormalon technique utilizes the fact that there is only one dimensional parameter in QCD, 
that is Aqcd- Models with dominance of a single chain introduce explicitly or tacitly other 
mass scales. In the cases mentioned above it is either the gluon mass or the scale where the 
coupling freezes out. 

3 In fact, as is argued in Ref. [0, two- loop effects are not suppressed within this model, while higher 
loops are negligible. For simplicity of presentation, we do not discuss this complication. 
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The phenomenology of the power corrections is rich enough to distinguish between various 
models. In particular, let us concentrate on the thrust T (see section 2) and heavy and light 
jet masses, M% andMf, respectively Then, the one-chain dominance implies || [|: 

<1-T> l/Q » (^>i/Q, (28) 
.M 2 .M 2 

\-Q2h/Q < \qj)i/q- 

Indeed, in the one loop approximation, one jet becomes "heavy", M\ ~ a s (Q 2 )Q 2 while the 
other jet remains "light", i.e. represented by a bare quark. Moreover, the pattern of the 
power-like corrections repeats the pattern of one-loop perturbative terms since it is a part 
of this correction. The latter statement i s especially obvious in case of the non-zero gluon 
mass. 

To summarize, the model with freezing of the running coupling, as applied to the power 
corrections, is in variance with the conclusion on importance of the multi-renomalon chains. 
Moreover, it leads to a well defined pattern of power corrections and can b e tested experi- 
mentally. We shall discuss this point further in Section 3. 

Ultraviolet renormalon. 

So far we considered contributions of momenta much smaller than the typical mass scale 
Q. Integration over virtual momenta brings contributions of momenta much larger than Q 
as well. These are detected through ultraviolet renormalons which are discussed in more 
detail in this subsection. 

Imagine that we are considering a Feynman graph which is convergent and saturated 
by typical momenta of order Q. If we are interested specifically in the contribution 5jjv of 
virtual momenta k 2 > Q 2 , then we could start with the following estimate: 

Q ^-*s{k 2 ) (29) 

where we put the factor Q 2 to make the quantity dimensionless. For example, we may always 
consider the fraction of the total answer which is due to very large momenta. Note, however. 
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that because because the coupling depends on k 2 only logarithmically, the integral (]29D is 
saturated at its lower bound by k 2 ^ ~ Q 2 and we would end up with consideration of the 
same typical momenta of order Q. 

The situation is changed if we substitute the expression for a s (k 2 ) as an expansion in 
a s (Q 2 ) and consider large-order terms in this expansion. Then for the coefficient a n of the 
perturbative expansion in a s {Q 2 ) we get the integral: 



which is dominated by k 2 jj ~ e n Q 2 , where e is the base of the natural logs. Thus in this 
case k 2 fj- Q 2 and we encounter the ultraviolet renormalon. 

To get a rough estimate of the effect of two renormalon chains we consider an expression 
similar to (|29|) but with a s (k 2 ) replaced by a 2 (k 2 ). Applying again Eq. (f26|) we can readily 
see that two renormalon chains result in the same asymptotical behaviour of the coefficients 
a n at large n as a single chain. 

Such calculations can be performed in much greater detail J7|. In particular two- loop 
effect can be accounted for systematically, also in the /3-function. The results are in agree- 
ment with the estimates above. Moreover, one can show that the "towers" of momenta, 
/c| 3> k\^> Q 2 still result, generally speaking, in the same asymptotical behaviour of a n as 
above. Instead of going into technical details, let us mention another simple way to under- 
stand the effect. Namely, one might argue that the renormalon chain is singled out because, 
to a given order n, the power of the log in the integrand (see Eq. (|3"0"D) is the highest possible. 
Moreover, since logs are eventually large, ln{k 2 jj/Q 2 ) ~ n this contribution appears to be 
the largest. Indeed, two renormalon chains produce one power of the log less. However, 
the two-chain contribution is enhanced combinatorially. The reason is that now the same 
number of the bubble insertions, n, can be distributed between two chains. As is readily 
seen, this combinatorial factor is also of order n. Thus, two mechanisms of the growth of a n 
with n, namely large combinatori al factors and large logs in fact get mixed up. 




(-1) 




(30) 
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Note that so far we have been considering the perturbative expansion in a s (Q 2 ) which is 
most natural. However, one might consider expansion in the coupling normalized at arbitrary 
point /i 2 as well. For example, expansion in the bare coupling a s {A 2 JV ) is relevant to the 
lattice studies (with huv ~ a -1 ) j25| - Then the perturbative series is reshuffled considerably. 
In particular, in terms of the expansion in a s (/i 2 ) the ultraviolet renormalons are power-like 



suppressed if /i 2 ^> Q 2 p6f and disappear completely if /i 2 = Afjy. Indeed, there is no 
integration over virtual momenta k 2 > Afjy at all. Instead of (p0|) we now have: 



Q 2 a s (k 2 )— « if- £ (&o) n ^<(A 2 a y ) ,iVg = ln{A 2 uv /Q 2 ). (31) 



Thus, we have a power suppressed, as Q 2 /A 2 JV , term with a series in a s (A^ v ) in front which 
behaves as an infrared renormalon in an unusual position and with n! behaviour tempered 
at n « Ag. Thus, in terms of the expansion in a s (A 2 /v/ ) the ultraviolet renormalon appears 
to be related to the power suppressed terms of order Q 2 /A 2 rv . There are first indication s on 
the existence of terms ~ (Aq C d • a ) 2 in the lattice simulations but the correspondence 
(if any) with UV renormalons has not been considered in any detail. 
Effective four-fermion interaction. 

The UV renormalon is most effectively treated within an operator product expansion 
which utilizes the observation k 2 ^ ^> Q 2 | f2~7j ^, |28 |. Moreover, since the ultraviolet renor- 
malon is related to very high momenta, the structure of this OPE is universal. The leading 
operators are singled out through their anomalous dimensions and various channel differ only 
in values of matrix elements of these operators . 

It is remarkable that at the level of two chains there appear four-fermionic operators 
which actually dominate the large n behaviour of the expansion coefficients a n [0, |29|, ^8 



For purpose of orientation, one may keep in mind that the dominant operator Oq is of the 
form: 

c(k 2 )0 6 = Ew ^ (fc2) (qt b qqt b q + qt b l5 qqt b l5 q) . (32) 
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where t b is the Gell-Mann matrices in flavor space and k 2 is the effective k 2 ^ associated with 
ultraviolet renormalons, k 2 ^ ^> Q 2 (see above). Moreover, c(k 2 ) is a coefficient function 
which can in principle be obtained as an expansion in a s (k 2 ). Although there are many sim- 
ilarities of the technique with the standard operator product expansion, there are important 
differences as well. In particular, the matrix elements of the operators can be calculated 
perturbatively as far as Q is large while truncation of the expansion for c(k 2 ) at any finite 
order in a s {k 2 ) changes the resulting coefficients a n by order unit (for details see 0). 

It is remarkable that similar four-fermion effective interactions are introduced in Nambu- 
Jona-Lasinio model, with k 2 replaced by a mass scale of order Aq CD . It is for the first time 
that a four-fermion effective interaction is generated by well defined graphs in fundamental 
QCD and one may speculate further on matching of the UV renormalon with the NJL model 
at an intermediate mass scale [f3C[| . 

Dispersive approach to the running coupling: 1/Q 2 terms. 

The basic feature of the ultraviolet renormalon is that the corresponding perturbative 
series is Borel summable. Indeed, a simple function of the running coupling: 



has the same perturbative expansion as produced by the ultraviolet renormalon. The func- 
tion (|D is the Borel sum of this expansion. Note that the sign alternation of the expansion 
coefficients, a n ~ (— l) n is crucial for the Borel summability. 

If a n ~ (— b ) n n\ then the absolute values of the subsequent terms in the perturbative 
sum fall down until a critical value N rr is reached: 



For n > N cr the values of \a n (a s (Q 2 )) n \ grow with n. The Borel summation replaces this, 
growing branch of the perturbative expansion by terms of the order aN cr {<y s {Q 2 )) Ncr ■ In other 




(33) 



1 



(34) 



N, 
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words, the Borel summation effectively introduces a nonperturbative term of order 1/Q 2 : 

A 2 

~Q 



^non-pert ~ — ' (^) 



Although such a term can be small numerically and subordinate to many terms in the 
perturbative expansion with n < N cr , it brings in a problem of interpretation ||. Indeed, 
in the spirit of the quark-hadron duality the 1/Q 2 correction would be naturally associated 
with contribution of low-lying hadronic states since the corresponding imaginary part is 
proportional to 5(Q 2 = 0) and should be smeared over Q 2 ~ Aq CD to imitate the hadronic 
contribution. On the other hand, as we saw in the preceding subsections the ultraviolet 
renormalon is associated with huge virtual momenta, k 2 ^ ^> Q 2 . Thus, the duality looks 
unnatural. Note that in case of the infrared renormalons the s ituation is very different and 
natural: IR renormalons are dual to low-lying hadronic states. Similar basic idea underlies 
the QCD sum rules [[H]] (for a review at this Workshop see Ref. |32|]). The ultraviolet 



renormalon defies this logic by introducing duality between low-lying hadronic states and 
very large virtual momenta in perturbative graphs. 

This kind of connection between the imaginary part at low Q 2 and real part at asymptot- 
ically large Q 2 , arises if one postulates dispersion relations. In the context of the ultraviolet 
renormalon, the dispersion representation is to be used for the running coupling itself. Then 
1/Q 2 terms are introduced in the definition of coupling [33, |19| (analytical properties of the 



running coupling can be of course discussed without any connection to renormalons, see, 
e.g., Refs. 

The connection between IR and UV momenta is realized then in the following way. In 
perturbation theory, the coupling satisfies dispersion relations and the corresponding cuts 
are at s > 0. If we assume analyticity for a resummed coupling as well then the standard 
expression, a s (Q 2 ) ~ 1 / b ln(Q 2 / Aq CD ) , contains a pole at Q 2 = Aq CD . On the other hand, 
no finite order of perturbation theory has a c ut below s = 0. The non-vanishing imaginary 
part at the tachyonic momenta is thus introduced via the summation procedure of all the 
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orders in a s (Q 2 ). The summation procedure is well known to be justified in the leading log 
approximation and is perfectly consistent at this level. 

It is an open question how to define the coupling to the accuracy Aqcd/Q 2 - In particular, 
one may introduce a running coupling with the Landau ghost removed 



2 



asiQ2) b ln(Q 2 /A 2 QOD ) + bo(A 2 Q cD- Q 2 ) ' (36) 



At large Q 2 the coupling a s (Q 2 ) differs from a s (Q 2 ) by a negative 1/Q 2 correction. Note 
that it is the tacit assumption on the analyticity that relates the removal of the Landau pole 
to the 1/Q 2 correction at large Q 2 . 

The crucial point which we are emphasizing again and again in this section is that higher 
loops in perturbative expansion collapse to the same order of magnitude power corrections. 
Within the current set up of the analyticity, consider the effect of removing the Landau pole 
from a™(Q 2 ). Note that we should remove the single pole: 

1 1 ^QCD (37) 



b%ln"(Q 2 /A 2 QCD ) b$ln"(Q 2 /A 2 QCD ) (n - 1)%{Q* - A 2 QCD ) ' 

and we come again to a 1/Q 2 correction which is not suppressed by any power of a s (Q 2 ). In 
particular, if we have an ultraviolet renormalon then 

£n!(-M n < +1 W 2 ) - -^E(-l)" (38) 

as far as the 1/Q 2 corrections are concerned. To make sense out of the latter sum we can 
introduce again the Borel summation (which sums it up to 1/2). In this way we see again that 
the Borel summation of the UV renormalon and introduction of the Aq CD /Q 2 corrections to 
the coupling are interrelated. A condition for the consistency of the two approaches might 
well exist but has not been investigated. 



Eq. (|37D demonstrates [[y| that the procedure of the removal of the Landau pole does 
not reduce to a universal redefinition of the coupling on the 1/Q 2 level (as is usually tacitly 



20 



assumed). Instead, the l/Q 2 correction depends on all the terms in perturbative expansions 
in a nontrivial way. 

Summary to part II. 

To summarize, we have demonstrated that the power suppressed terms receive contribu- 
tions of the same order from higher loop graphs. The observation can be phrased different 
way as well. Namely, in case of the logarithmic terms (pure perturbative QCD) larg e logs 
are known to cancel between virtual and real gluons as far as inclusive observables are con- 
cerned. Because of these cancellations higher loops are suppressed by extra powers of a small 
coupling a s (Q 2 ). No such cancellations are known for the power- suppressed corrections. 

The accumulation of the contributions due to higher loops could provide a pure pertur- 
bative mechanism for an enhancement of the power corrections. However, it is not known 
at all whether this enhancement is actually taking place. And even if one is prepared to 
speculate that the enhancement is there, experimental consequences of this hypothesis are 
not clear since there is no model independent way to relate various channel. We will consider 
some models in the next section. 



4 Dynamical models. 

Infrared renormalons and hadronization models. 

It seems rather obvious that the IR renormalons should converge with simple hadroniza- 
tion models. Indeed, if we look afresh at, say, Eq. fllQ|) then we immediately recognize 
renormalons as introduction of "intrinsic" transverse momentum for partons, or the tube 
model (for a review see Ref. [PH]). Moreover, the same l/Q corrections were first discussed 



within this kind of an approach long ago | 36fl . However, it is not straightforward at all to 



quantify this feeling and fix the relation between renormalons and hadronization models in 
any rigorous way. The problem is the same as in the previous section: on one hand, all 
higher loops become equally important for the power corrections and, on the other hand, 
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there is no han die to treat all the terms in perturbative expansions. 

The statement on the universality of the 1/Q corrections || [nj falls closest to fixation of 



the relation between the renormalons and hadronization models. To derive the universality 
one ensures first that nearly two jet configurations dominate an observable. This can be 
achieved in a p ure kinematical way. For example, in case of the thrust one introduces the 



following average [[L0 



(exp{-u(l-T))) 1/Q = exp{-2vU/Q) (39) 

where v is a large number and U parameterizes infrared region. The dynamical question is 
how to find U. 

In the universality framework and to first order in a s (k\): 

u = %r Qd ^* s {ki) ~ (40) 



-kQ Jo k± Q 



Although Eq. (40) looks very similar to the single chain approximation they are not identical. 
Namely, there is a crucial factor of 2 in the value of U. The justification for this factor of 2 is, 
roughly speaking, as follows. We start with perturbative graphs and pick up the contributions 
which contain logs and dominate the thrust distribution (see the remark following Eq. (|5|)). 
Then we continue these contributions to the infrared sensitive region using the renormalon 
technique. The factor of 2 is due to the fact that originally there are large contributions 
associated with gluons aligned with each quark. Note that the momentum configuration 
which finally determines the 1/Q correction corresponds to a soft gluon, not a collinear one. 
Hence, there is no solid proof of this factor of 2 |J7| . 



With inclusion of higher loops, the universal factor U becomes 0, 10 



r~ Q dk 2 

U = / y^JeikMkl)), (41) 
Jo Qk± 

where 7^ is the so called cusp anomalous dimension. This substitution allows for extrap- 
olation to the infrared- sensitive region of all the log terms in the perturbative calculations. 
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Since these logs are known to be universal the power corrections turn universal as well. The 
universality of the 1/Q corrections implies, in particular, that the same quantity U enters 
the expressions for other shape variables, e.g. ||: 

\{l-T)y Q = ^(C) l/Q = -(-)i/Q = -{EsinH) 1/Q = U. (42) 

Here we used the standard notation for various shape variables and for simplicity did not 
introduce the exponential weighting (see Eq. (|39|)). 

From the phenomenological point of view the predictions from the universality of the 
1/Q corrections are close to the tube model. In particular, we expect 0: 

,M?, ,M 2 . , 

(l-r>l/Q = (gf)l/Q+(g^)l/Q (43) 

Ml Mf 

which is to be contrasted with an alternative set of predictions (28) based on the assumption 
on the dominance of a single renormalon chain. It is our understanding that the data |38| 
clearly favor the multi-renormalon, or universality models but the fits seem to be subject to 
large systematic uncertainties. 

UV renormalons and short distance physics. 

While everybody is prepared to consider models for infrared phenomena, treatment of 
the ultraviolet region is expected to be rigorous. Indeed, the effective coupling is small 
and the nonperturbative effects arguably vanish as a high power of Q^ 1 Ultraviolet 
renormalons are associated, on the other hand, with Aq CD /Q 2 corrections. Thus, if we 
assume that UV renormalons give indeed a significant contribution then we introduce in fact 
new nonperturbative effects at short distances. Let us first substantiate this simple point. 

Consider to this end again the heavy quark potential at short distances. If we allow for 

a 1/Q 2 piece in the effective coupling (see Eq. ( Pop above) then we have for deviations of 

the potential from the Coulombic one at short distances: 

, a 2 

lim 5V(r) ~ d 3 kexp(ik-r)^^ ~ A 2 QCD ■ r. (44) 

T ► J j£ 
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Such a linear piece does not arise from the infrared region at all. Moreover, the heavy quark 
potential is measurable on the lattice (see, e.g., |EJ) and in this sense is observable. Let us 
emphasize that we are discussing short distances while the linear potential is of course very 
common at large distances but is expected to change to a r 2 correction at short distances 
(see Eq.(p3|)) according to the conventional picture (see discussion in Section 3). So far, this 
prediction has not been checked. Moreover, there are first indications of the 1/Q 2 correction 



to the effective coupling obtained on the lattice pqj . Thus, it might be right time to look 
for implications of the 1/Q 2 corrections associated with short distances. 

To visualize what kind of new physics we can deal with let us consider a simple estimate of 
nonperturbative effects in V(r) at short distances Namely, let us represent the potential 



energy of the QQ pairin an abelian case as an integral over space from the quark electric 
field: 

V( r ) ~ -L />r'Ei(r') .E 2 (r + r'). (45) 

47T J 

As far as we substitute the electric fields E 12 as the fields of point-like charges, the Coulomb 
potential is perfectly consistent with (|45| ) of course. Imagine, however, that the electric field 
is changed at large distances because of the confining properties of the vacuum. Then there 
would be a corresponding change in the potential at short distances as well. Moreover, let 
us assume that strong perturbative fields are allowed in the vacuum while weak fields below 
some critical value are not allowed: 

o 

cy r 

Kr ~ Kcd or -j- < Aq CD (46) 

^cr 

where we accounted for the fact that at the distances large compared to r the field of the 
QQ pair is that of a dipole. As the last step, estimate the change in the potential at short 
distances due to the screening of the field at r > R cr : 

im^m ~ ^Cw ~ a i ~ al/2Ahc ° r ' (47) 

and we reproduce the linear correction to the potential at short distances. 
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It is worth emphasizing, however, that R cr ~ r 1 / 3 and tends to zero with r — > 0. Thus, 
the dynamics described above implies that the vacuum is fine grained. Which is far from 
being obvious. Alternatively, we may say that if the linear correction to the potential at 
short distances is confirmed this would have highly nontrivial implications for QCD. 

Tachyonic gluon mass. 

If we decide to look for a dynamical framework for a new short- distance physics in QCD, 
then monopoles seem to be a natural candidate (the topic of monopoles was extensively 
covered at this Workshop and for reviews see, e.g. The monopole s appear in the 

abelian projection and are usually discussed in connection with the confinement mechanism, 
or physics at large distances of order Aqq D (see, however, pE2]|). Here, we are interested in 
the short-distance aspect of the problem. Moreover, we would argue that if the monopoles 
are indeed an important field configuration in QCD, then a new short- distance physics is 
introduced. Indeed, in theories with elementary Higgs fields, the singular behaviour of gauge 
fields at short distances is smeared. However, in QCD at short distances no composite field 
can be relevant and the singular behaviour of the gauge fields can be regularized only at 
the lattice size, at least at first sight so. These are very interesting questions to our mind 
but they were not discussed, to our knowledge, in original papers and we will, therefore, 
confine ourselves to a few remarks sketching a possible new kind of phenomenology with a n 
effective tachyonic gluon mass, —A 2 < 0. This time the effect is coming from short distances 
and the A 2 is not to be confused with A 2 considered in Section 2 in connection with the IR 
renormalon (or non-analytical in A 2 terms). However, there is a common feature of A 2 and 
A 2 : in full nonabelian theory both can be consistently tried only in one- loop calculations. 

The abelian models with condensation of monopoles can be viewed as (dual) Higgs mod- 
els, with two masses m 2 H and A 2 , of spin and spin 1 particles, respectively. The quark static 
potential is then given by (see, e.g., 



Q 2 e~ Xr Q 2 A 2 /A 2 + m 2 \ 

V(r) = - -f + In rr-^ • r (4£ 

v ; 4vr r 8vr \ A 2 ) v 
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where Q is the abelian charge. Moreover, the first term represents Yukawa-type potential due 
to a (massive) gluon exchange while the second term represents the confining potential kr, 
where k is the string tension. As is mentioned above, the potential Q38| ) is considered usually 
at large r. For us, the central question is how far we can extend it to the region of small r. 
The answer is easy to obtain in the the London limit, m 2 H ^> A 2 . Then the linear piece in 
( fi8|) can be trusted also at A -1 < r < m^ 1 . This can be seen from more detailed explicit 
expressions for the potential or one may simply observe that in the limit mu — > oo 
the core of the string connecting the quarks becomes infinitely thin. Let us mention that 
there is numerical evidence jPJ] that the realistic QCD is close to the London limi but this 
is far from being a settled question (for further discussion see j4"Ifl). We simply assume 
mil ~ A[7y since it is the easiest way to implement zero siz e of monopoles (see discussion 
at the beginning of the subsection), although other possibilities may also exist. 
Then for the correction to the potential at short distances we get: 

]im5V(r) = -X 2 r(^-^ln(^-\] = - t^r. (49) 
r- o v ; \8ir 16vr \ A 2 J J 2 K J 

Here A 2 is positive and the expansion corresponds to an effective tachyonic gluon mass —A 2 . 
Numerically, A 2 ~ 3k/2a s . Note that our simple-minded estimate ([|7]) also corresponds to 
a positive 5V(r) and, therefore, could be interpreted as due to a tachyonic gluon mass. 

We could have tried to scrutinize the derivation but would prefer to jump directly at the 
phenomelogical manifestations of the tachyonic gluon mass. In fact the phenomenology turns 
not trivial at all and, to our knowledge, does not echo any other approach . We illustrate 
this statement by example of the QCD sum rules in the p-meson and pion channels. As is 
well known ||5] these channels exhibit striking differences and, essentially, the pion channel 
cannot be understood in the standard terms of the IR sensitive power corrections. 

If one accounts for a non- vanishing gluon mass A then in one-loop approximation the ab- 
sorptive parts of the current correlators in the vector and pseudoscalar channels for massless 
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quarks look as 0: 

ImUJs) = 1 + — (50) 

7T 

. v 17a s 2a s , u 2 Aa s A 2 

imIL, s = 1 + — - + — In" - — . 

67V 71 S 7T S 

Here \x is the normalization point and the term proportional to Zn(/i 2 /s) is due to a non- 
vanishing anomalous dimension of the pseudoscalar current. Moreover, the currents are 
normalized in such a way that for bare quarks (a s = 0) the imaginary part is equal to unit. 
Note also that the gluon condensate would manifest itself as a term proportional to S' (s) 
while in the equations above s ^ 0. 

We are interested now in the effect of the gluon mass which affects only the pseudoscalar 
channel and in accordance with our discussion above we substitute a tachyonic mass instead 
of A 2 so that a rough estimate runs as follows: 

4a s A 2 2k OAGeV 2 
► « + . (51) 

7TS S S 

Thus, the effect reaches ~20% at s ~ 2GeV 2 which fits very well, both in the magnitude and 
sign, what is needed to bring the sum rules in agreement with the data, the agreement which 
is out of reach |^| of the conventional phenomenology f\ It is amusing that the correction 
~ A 2 does not appear in the vector channel where there exist quite strict bounds on the l/s 
terms |46|| . Preliminary analysis shows that the A 2 could resolve the problems |^5| with sum 
rules in the gluonic channels as well. It is worth mentioning that the preliminary evidence 
for the 1/Q 2 correction on the lattice []25| could also be fitted with a tachyonic mass for the 
gluon. 

Summary to part III. 

In this part we have considered dynamical models for nonperturbative effects in QCD, 

both in infrared and ultraviolet, motivated by multi-renormalon chains. In case of the IR 
4 I am deeply indebted to K.G. Chetyrkin for providing me with these equations. 

5 There exist successful fits within the model of the instanton liquid E7|] (see also |45[|). It is not neces- 
sarily conflicting models since monopoles may be interconnected with the instantons (see [}4l| and references 
therein). 
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renormalons we argued in favor of convergence of the renormalons with the old fa shioned 
tube model. In case of the UV physics, we speculated about applicability of the string 
model with infinitely thin core. Phenomenologically, the model amounts to introduction of 
a tachyonic effective mass for gluons and there are first indications that the model might 
resolve old standing problems with QCD sum rules. 
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